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Abstract 

Density functional theory and molecular dynamics simulations have been used to optimize 
the structure of nanowires of Si02- The starting structures were based on j3-cristobalite, ortho- 
tridymite, j6-tridymite, and rutile crystals. The analysis of the electronic structure has been 
validated by many-body perturbation calculations using the Go Wo and GW + Bethe-Salpeter 
equation approximations, in order to account for quasi-particle and excitonic effects. The 
calculations indicate that many of these nanowires have semiconducting character, while the 
corresponding bulk solids are insulators. In the case of thick rutile-like nanowires we found 
the gap congested with surface states. Electronic charge is transferred from silicon atoms 
to oxygen atoms, and the bonding in the nanowires is partly ionic and partly covalent. The 
magnitude of the band gap can be engineered by changing the structure and the thickness of 
the nanowires. 

Keywords: density functional theory; silicon dioxide; nanocables; semiconductor; surface states; 
many-body corrections. 



INTRODUCTION 

Metallic oxide uses, including as catalysts, and in state of the art batteries, fuel cells, solar cells, 
and supercapacitors for electronics^^ is recently arising as a burgeoning area, especially for 
nanostructures, like nanoclusters, nanotubes, nanowires and nanocables Recently, it has been 
shown that a metallic oxide combination, Ru02/Si02, consisting of a transparent conductive plat- 
inum group metal oxide and a dielectric insulating silica material, can be made in a cylindrical lay- 
ered fashion J-^ The Ru02 part forms a thin shell around the inner solid core nanowire (or nanorod) 
of Si02- Because Ruthenia is such an extraordinary material for charge storage - the highest known 
- being able to use cores of an inexpensive material to act as its scaffolding becomes an ideal way 
to use the expensive metal ruthenium. 

Work on determining theoretically the intrinsic quantum conductances and capacitances,^ as 
well as the junction capacitances,^ of ruthenium oxide nanowires and nanocables, followed the 
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earlier effort in producing conductive Ruthenia nanoskins on insulating silica paper.^The essen- 
tial purpose of those studies was to investigate charge storage for chemical and electronic appli- 
cations. What must be done now is to completely characterize the underlying Si02 nanowires, 
since their properties may affect the overlying RuC>2 material. Thick SiC>2 nanowires have been 
synthesized.^^ However, to our knowledge, few theoretical studies of SiC>2 nanowires have been 
performed, and only for very thin nanowires J 17 * 18 ! We here present density functional calculations 
of the atomic and electronic structure of SiC>2 nanowires as a function of nanowire thickness. An 
analysis of the electronic density of states, electrical character, electron density, charge transfer 
and chemical bonding is presented. 
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Figure 1: (Color online) Two views of the crystalline structure of some Si02 bulk-phases: A) /3- 
cristobalite;^ B) ortho-tridymite; 30 C) /3-tridymite.^ Oxygen and silicon atoms are represented 
by red and tan spheres, respectively. 



METHOD 

In the density functional theory (DFT) calculations of the structure of the Si02 nanowires we have 
used two different simulation packages: the localized basis set code Fireball, 19 used for the 
full-characterization of the structures, and the plane-wave code D AC APoJ 20 * 21 1 used for refining the 
electronic properties of the optimized structures. The computational Fireball scheme, as well as 
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its theoretical foundations, has been described in full detail elsewhere. 

| 19 | 22 | 23 | FIREBALL 

is a super- 
cell implementation of DFT formulated in real space using a basis of local pseudoatomic-orbitals. 
Here we only summarize the main points. A self-consistent version of the Harris-Foulkes func- 
tio nal™ is used for solving the Kohn-Sham (KS) equations where the KS potential is calculated 
by approximating the total charge by a superposition of spherical charges around each atom. The 
local density approximation (LDA) is used for the description of exchange and correlation (XC) 
effects. An extended sp 3 d 5 basis set of single numerical atomic orbitals (one s, three p and five d 
orbitals) is employed for Si, and a double numerical sp 3 s*p* 3 basis set for O. Core-electrons are 
replaced by norm-conserving scalar relativistic pseudopotentials.^ 

In the DACAPO code,' 2 ^^ the KS equations are solved using a periodic supercell geometry. 
The super-cells are such that the separation between neighbor nanowires is at least 15 A, in order 
to avoid inter-nanowire interactions. The ion-electron interaction is modeled by ultra-soft pseu- 
dopotentials,^ and XC effects are treated by the generalized gradient RPBE functional. 20 The 
one-electron valence wave-functions are expanded in a basis of plane- waves, with energy cut-offs 
of 400 and 500 eV for the kinetic energy and the electronic density, respectively. The supercell 
size and the energy cut-offs have been adjusted to achieve sufficient accuracy in the total energy. 

The Brillouin zones of all the bulk-phases and nanowires were sampled with 6x6x6 and 
12 x 1 x 1 Monkhorst-Pack grids ^respectively. In order to obtain the minimum energy structures 
with Fireball, a dynamical quenching technique has been employ ed.^^ Within this approach, 
all atoms are allowed to relax until an energy minimum is reached. Then, the atoms are displaced 
slightly and the dynamical quenching process is continued. We run the dynamical relaxation pro- 
cess until the force acting on each atom was less than 0.01 eVA~ 1 . At the same time, we minimized 
the strain on the unit cell in all periodically repeated directions. 

Once the optimized structures have been established with the Fireball code, these structures 
are used as starting geometries in a steepest-descent re-optimization, carried out with the DACAPO 
code. For all the systems, the convergence of the optimization procedure was immediate, revealing 
the consistency between the two optimization approaches. Additionally, in order to check the 
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Figure 2: (Color online) Front and lateral views of the most stable Si02 nanowires obtained by 
relaxing fragments of bulk-phases (see [I]): A) /3-cristobalite; B) ortho-tridymite; C) /3-tridymite. 
Front views of the initial nanowires before relaxation are also shown. Oxygen and silicon atoms 
are represented by red and tan spheres, respectively. The black dashed boxes indicate the unit cell 
along the nanowire axis. 
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Figure 3: (Color online) Lateral and front views of the most stable Si02 nanowires with rutile- 
like structure. Views of the initial nanowires A and B before relaxation are also shown (Ao and 
Bo). From top to bottom (A to F) the nanowires have 4, 12, 16, 24, 32 and 40 Si02 functional 
units per unit cell, respectively. Oxygen and silicon atoms are represented by red and tan spheres, 
respectively. The blue dashed boxes indicate the unit cell along the nanowire axis. 
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Figure 4: (Color online) Cohesive energy of nanowires per Si02 functional unit (as defined in the 
main text) as a function of the number of functional units in the unit cell, for the nanowires shown 
in [2] and [3) 
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influence of the temperature on the structure of the nanowires, molecular dynamics simulations at 
temperature T=500 K have been performed on some of the most stable structures, which revealed 
the high structural stability of the nanowires. In all cases, a typical dynamical pattern of two 
vibrational modes 10 is observed: a radial, breathing mode and a longitudinal, stretching mode, 
conserving the total energy in time (in tests up to 10000 ps). 

GEOMETRICAL STRUCTURE OF THE NANOWIRES 

Due to the lack of information in the literature about the possible atomic structures of small di- 
ameter SiC>2 nanowires we have approached the problem from a heuristic point of view. As a 
reasonable initial choice we constructed SiC>2 nanowires based on one dimensional-like fragments 
of the most stable bulk phases. A number of polymorphs of silica, SiC>2, exist in nature. Co- 
esite and stishovite are polymorphs which form at high pressure and temperature. At moderately 
lower temperatures one finds the tridymite and cristobalite crystalline forms. These, like quartz, 
present low and high-temperature structural modifications. The crystalline lattices of three of those 
phases, /3 -cristobalite, 29 ortho-tridymite, 30 and /3 -tridymite^ are plotted in[TJ According to this, 
we have adopted as starting nanowire geometries fragments of the /3 -cristobalite, ortho-tridymite, 
and /3-tridymite crystals in the form of SiC>2 nanowires. Those fragments are of infinite length 
in the direction of the nanowire axis. The unit cell contains six SiC>2 units in the case of the /3- 
cristobalite nanowire, and twelve SiC>2 units in the two tridymite wires. Dynamical relaxations led 
to the optimized structures of [2j where frontal and lateral views of the most stable nanowires are 
shown. [2] also shows a front view of the starting bulk- like nanowires before relaxation. Molecular 
dynamics simulations revealed that the minimum energy structures obtained are robust and stable 
up to 500K at least. Optimization of the cell parameter in the axial direction was also allowed, 
but no significant variations were observed with respect to the starting cell parameters. The cal- 
culated lattice parameters are 5.07 A, 8.27 A and 8.24 A for the /3 -cristobalite, /3-tridymite and 
ortho-tridymite-based Si02 nanowires, respectively. The starting bulk-like nanowires present a 
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hexagonal-like section with some oxygen atoms popping out due to stoichiometry. The section 
deforms after structural relaxation: the hexagonal inner cavity of the nanowire reduces its size and 
more oxygen atoms appear on the outer part of the nanowire. The reconstruction is more severe 
in the two tridymite nanowires. We also did calculations in which the wires were cut from the 
crystal choosing other directions for the wire axis. However, those led to less stable structures 
after optimization. 

Table 1 : Average nanowire diameter estimated from the most external radially opposite Si atoms, 
Dsi, average diameter estimated from the most external radially opposite O atoms, Do, and average 
Si-0 bond-length, dsi-o, m the bulk-based and rutile-based Si0 2 nanowires. All distances are 
given in A. 





D Si (A) 


D (A) 


dsi-o (A) 


Bulk-based 


/3-cristobalite 


5.03 


6.73 


1.73±0.06 


Ortho-tridymite 


6.31 


8.74 


1.74±0.08 


/3 -tridymite 


6.43 


8.45 


1.74±0.08 


Rutile-based 


4 Si0 2 units 


4.67 


6.57 


1.73±0.02 


12 Si0 2 units 


10.32 


12.23 


1.74±0.05 


16 Si0 2 units 


14.12 


15.64 


1.74±0.04 


24 Si0 2 units 


16.09 


17.84 


1.73±0.05 


32 Si0 2 units 


18.67 


21.16 


1.74±0.06 


40 Si0 2 units 


23.04 


24.88 


1.74±0.08 



We have performed additional calculations starting with other structures. For this purpose we 
first took a narrow Si0 2 nanowire whose lateral and front views are shown in panel Ao of [3] That 
structure, having four Si0 2 molecules per unit cell, is taken from previous works's^ on nanowires 
of transition metal dioxides (M0 2 ), and it was obtained by relaxing nanowires cut from the rutile 
crystals. Shown in panel A is the relaxed Si0 2 nanowire, whose structure is more symmetric than 
that of the starting nanowire. Next we relaxed a nanowire with the form of a two-wall nanotube, 
whose lateral and front views are shown in panel Bo of [3] This structure was also taken from 
previous work on transition metal dioxide nanowires.^ Structural relaxation distorts the original 
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nanotube producing an atomic arrangement based on the rutile lattice (panel B). Motivated by 
this feature, we explore a growth pattern to construct thicker nanowires based on the rutile-like 
structure of panel A. Nanowires with 16, 24, 32 and 40 Si02 units per cell are shown in panels C, 
D, E and F. These structures are shown after relaxation. For nanowires of a given size, different 
starting configurations were considered, but only those leading to the most stable nanowires are 
shown in [3} The structures of all these nanowires are stable in molecular dynamics simulations at 
500K. 

The cohesive energies of these nanowires per Si02 unit, defined as: 

have been plotted versus the number of SiC>2 molecules per unit cell in|4j In (??), £'°'[NW] is 
the total energy of the nanowire per unit cell, N is the number of SiC>2 molecules per unit cell, and 
E tot [SiC>2 unit] is the total energy of an isolated gas-phase SiC>2 molecule. The cohesive energies 
of the reconstructed ortho-tridymite and /3-tridymite nanowires are similar, and higher than the 
cohesive energy of the reconstructed /3-cristobalite nanowires. That is, the tridymite nanowires are 
more stable. 

In order to check the validity of the structures obtained in the simulations using the canoni- 
cal (minimal) unit cell, we have carried out calculations taking a cell doubled in size, along the 
nanowire axis, compared to the original one. The aim of the test was to allow for the possibil- 
ity of observing new structural reconstructions not seen in the dynamical quenches based on the 
minimum unit cell. The test was made for all the bulk-based nanowires and four rutile wires. 
The results obtained do not reveal any significant variation of the structures with respect to the 
nanowires obtained by just using the minimum unit cell. 

The cohesive energies of the rutile-based nanowires are also included in |4} A similar value of 
the cohesive energy per SiC>2 unit is obtained for all the rutile-based nanowires, although it tends to 
decrease as the cross section of the nanowire increases. That is, the cohesive energy per Si02 unit 
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tends to decrease for thick nanowires. We think this is due to the fact that no stable bulk phases 
of SiC>2 with structures related to rutile exist at normal pressures. The only bulk phase related to 
rutile is stishovite,^ a high pressure phase with a mass density of 4.29 g cm~ 3 , much larger than 
the density of a-quartz, 2.65 g cm~ 3 . Consequently, thick rutile-like nanowires are not expected to 
form at normal pressure. Thin nanowires, like those studied here, present fewer problems because 
the structural constraints are weaker. Still, one can notice the enhanced stability of the rutile-like 
nanowires with respect to those obtained from the /3-cristobalite, ortho-tridymite and /3-tridymite 
phases. The gain in cohesive energy with respect to the tridymite-like nanowires with a similar 
number of Si02 units ranges between 1 eV and 2 eV per Si02 unit, and from 2.5 and 3 eV per 
Si02 unit with respect to the /3-cristobalite-like nanowire. 

To interpret the behaviour of the cohesive energies in [4] one has to notice that the range of 
cross sections for these nanowires corresponds to the non-scalable regime, well known in clus- 
ters,' 3 - 3 -' where the properties vary in a non-smooth way, often showing oscillations, as the cluster 
size increases. The present calculations suggest that the non-scalable regime occurs also in thin 
nanowires, and is responsible for the oscillations of the cohesive energy of the rutile-like wires, 
and for the initial decrease of the cohesive energies of the /3-cristobalite, ortho-tridymite and /3- 
tridymite wires. 

Basic structural measurements of a nanowire can include its diameter as well as its length. 
These parameters may be related to laboratory preparation of nanowires.^^ At the atomistic 
level, a characteristic parameter would be an average bond-length. Accordingly, [T] provides the 
calculated average diameter estimated from the most external radially opposite Si atoms, the 
average diameter estimated from the most external radially opposite O atoms, Do, and the average 
Si-0 bond-distance, dsi-o> in ai l tne bulk-based and rutile-based Si02 nanowires. Notice that the 
largest diameter obtained is that for Do in the case of 40 functional units per unit cell. Its value is 
about 25 A=2.5 nm, which, from the point of view of experimental synthesis, corresponds to a thin 
nanowire. Although Si02 nanowires used as a core in composite systems such as RuCVSiC^ have 
tipically been one order of magnitude thicker,^ there is no reason why smaller core nanowires 
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could not be processed and employed. Since the thicknesses of the nanowires studied here are 
on the non- scalable regime, where the properties may vary in a non smooth way, it is not easy to 
extrapolate the results of |4]to very thick wires. 



ELECTRONIC STRUCTURE 



Density of states (DOS) 
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Figure 5: (Color online) Density of electronic states (DOS) versus energy for the Si02 nanowires 
of [2] The Fermi energy £f is taken as reference. Upper, middle and lower panels correspond 
to nanowires based on the /3-cristobalite, ortho-tridymite and /3-tridymite crystals, respectively. 
Black curves give the total DOS, and red and blue curves represent DOS projected, respectively, 
onto the p states of Si and O atoms. Lateral views of the nanowires are included as insets. Oxygen 
and silicon atoms are represented by red and tan spheres, respectively. 

We begin our analysis by considering the nanowires obtained from the cristobalite and tridymite 
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Figure 6: (Color online) Density of electronic states (DOS) versus energy for the rutile-based Si02 
nanowires of [3} The Fermi energy £p is taken as reference. The number of Si02 units per cell and 
the average radius of each nanowire are indicated in the inset. 
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Figure 7: (Color online) Band- structure diagrams (referred to the Fermi energy level) for ru tile- 
based Si02 nanowires (see [3]) with 4 (panel A) and 24 (panel B) functional units along the most 
relevant symmetry line from r — > Z (along nanowire axis). 
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bulk-phases (see [2]). In addition to its intrinsic interest, the study provides insight on how struc- 
tural changes in the nanowires may affect their electronic properties, bonding character and charge 
transfer between atoms. Such information could prove useful when trying to form Si02 semicon- 
ducting nano structured materials, as well as when forming composite nanocables with overlying 
metal or semiconducting materials. 

[5] shows the total electronic DOS profiles versus energy for the three nanowires of[2j In order 
to mimic the experimental resolution,^ a Lorentzian broadening of 0.05 eV has been applied to 
the energy levels. The magnitude of the electronic gaps, measured as the difference between the 
energies of the highest occupied (HOMO) and lowest unoccupied (LUMO) molecular orbitals, 
is 1.69, 1.57 and 1.06 eV for the /3-cristobalite (top panel), ortho-tridymite (middle panel) and 
/3-tridymite (bottom panel) nanowires, respectively, revealing the potential semiconducting-like 
character of these nanostructures, especially if appropriate donor or acceptor atoms were available 
to reside in shallow states below the band edges within the bandgap. The calculated band-gaps are 
smaller than the gaps of the corresponding bulk crystals. The gap of bulk /3-cristobalite obtained by 
DFT is in the range 5.5-5.7 eV for LDA calculations, and increases to 10.1-10.3 eV if quasiparticle 
corrections are taken into account.^ A gap of 6.6 eV is obtained in gradient-corrected BLYP 
calculations.^ Our RPBE calculations give a gap of 7 eV. Turning to the tridymite crystals, the 
BLYP gap of /3-tridymite is 6.6 eV,^ and our RPBE gaps for ortho- and /3-tridymite are 5.8 and 
7.7 eV, respectively. That is, the bulk phases are insulators. A substantial lowering of the gap of the 
nanowires with respect to the bulk crystals is then obtained. The band-gap values of the nanowires 
are similar to those for some elemental bulk semiconductors like Si, or III-V binary and ternary 
semiconductor compounds, like GaAs, InGaP, and InGaAs, with varying stoichiometry.^22] 

Regarding the DOS profiles of the three nanowires we observe different band structure behav- 
iors: i) the /3-cristobalite nanowire (top panel) acts as a pure acceptor semiconductor-like material, 
able to give up mobile holes whose volumetric density is 

p = Uve -(zF-£v)/k B T ^ jj^ (2) 
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Figure 8: (Color online) Density of electronic states (DOS) versus energy for the /3-cristobalite 
(A), ortho-tridymite (B) /3-tridymite (C) (see [2]), and the smallest rutile-based (D) (see|3| Si02 
nanowires. The Fermi energy £p is taken as reference. In each panel, DOS profiles are shown for 
the structures before and after structural relaxation. 
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since £ v ~ £p (here e v is the top of the valence band, kg is the Boltzmann constant and U v is 
a numerical constant depending on the temperature T); ii) the ortho-tridymite nanowire (middle 
panel) acts as a potential donor semiconductor-like material, able to create mobile electrons with 
density 

n = Uce -^c-£ F )/k B T w jj^ (3) 

since £ c ~ £/r (here £ c is the bottom of the conduction band and U c is a numerical constant depend- 
ing on T); and iii) the /3-tridymite nanowire (bottom panel) acts as an intrinsic semiconductor with 
tiny carrier densities 

n = p = ^/lUTc e -^- £ ^/ 2kBT , (4) 
with the Fermi level located in the middle of the gap 

£ v + £ c + k B T In (U v /U c ) £ v + £ c 
£f = ^ W ~T~ (5) 

for r or U v /U c tt\, and with no preference to donate or accept electronic charge 

In [5] we also show the projected DOS onto the p-states of Si and O atoms. In the energy range 
plotted, the occupied part of the DOS is dominated by the p-states of O. This means that the oxygen 
atoms have accepted some electronic charge from the silicon atoms (details are provided below). 
The sum of the Si and O p-orbital contributions falls slightly below the total DOS, the difference 
being the smaller s-orbital contributions. 

Table 2: Transport gaps (in eV), calculated within standard DFT, quasi-particle GoWo and GW + 
BSE formalisms, for the rutile-like nanowires with 4, 12 and 16 Si02 functional units. 



Transport gaps (eV) 


Standard DFT 


quasi-particle GqWq 


GW + BSE 


Rutile-like nanowire 


4 Si0 2 units 


4.53 


4.80 


5.98 


12 Si0 2 units 


2.84 


2.98 


3.58 


16 Si02 units 


1.62 


1.73 


2.01 
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Figure 9: (Color online) Average charge transfer per Si atom, in e~, from Si to O atoms, in the 
Si02 nanowires of [2] and [3j 
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[6] shows the total DOS profiles versus energy for the most stable rutile-based nanowires of|3j 
As observed in the top panel of [6} thin rutile-like nanowires (with 4, 12 and 16 SiC>2 molecules per 
unit cell) show a semiconducting behavior, and the electronic gap is progressively reduced as the 
number of SiC>2 functional units in the cell increases. For thicker nanowires, those in the bottom 
panel, new peaks appear in the DOS near the Fermi energy, and this region becomes congested. 
The structure of the electronic levels in that region can be seen in |7J That figure shows the band 
structure of two nanowires, with 4 and 24 Si02 functional units in the cell, respectively. The left- 
side panel shows clearly the gap existing in the thin nanowire. On the other hand, the band structure 
of the thick nanowire (notice the different energy scale) shows that a bunch of states occupy the 
region of the gap. These can be interpreted as surface-like states. 

A comparison between [5] and [6] shows differences between the electronic gaps corresponding 
to nanowires of similar thickness but different structures. The /3-cristobalite nanowire (with four 
Si02 units per cell) has a gap of 1 .69 eV, while the gap of the thinnest rutile-like nanowire (six 
Si02 units per cell) is 4.93 eV. For nanowires with 12 Si02 units per cell, the differences are again 
evident, although not so large: the rutile-like nanowire has a gap of 2.91 eV, and the two tridymite 
nanowires show gaps of 1.57 and 1.06 eV. The conclusion is that by varying the thickness and the 
structure of the nanowires, these can be engineered to produce a desired value of the electronic 

gap- 
There is a reduction of the band-gap of the nanowires with respect to the Si02 bulk crystals. We 
ascribe the difference to the particular structure of the nanowires, with many unsaturated oxygen 
atoms on the surface. This is required by stoichiometry, two O atoms per Si atom, and it becomes 
reflected, as shown in[5J in the dominance of the occupied DOS near the Fermi level by the p-like 
states of oxygen. The reduced gaps also indicate that the nanowires will be more reactive than the 
usual bulk Si02 phases. 

The structural relaxation of the nanowires shown in [2] and [3] is well explained by looking at 
the DOS. [8] shows the DOS profiles of four nanowires in the original unrelaxed structure and after 
structural optimization. In the cases of /3-cristobalite (panel A), ortho-tridymite (panel B) and j8- 
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tridymite (panel C), the density of states in the unrelaxed structure shows several high peaks in 
the region near the Fermi energy. Those peaks correspond to electronic states located inside a gap 
of width about 5 eV, which reminds the gap in the bulk crystals. The high DOS near the Fermi 
energy indicates that the unrelaxed structures of the nanowires are unstable. In fact, as the initial 
structures of [2] relax to more stable arrangements, the states at the Fermi level almost disappear, 
and a true gap develops, although smaller in magnitude than the original gap of the crystal. At 
difference with the previous wires, the smallest rutile -based nanowire shows a gap and no peaks at 
the Fermi energy. The effect of structural relaxation is, in this case, to widen the gap by moving 
the unoccupied electronic states up in energy. The widening of the gap increases the stability of 
the nanowire. However, as discussed above, surface states appear in the gap for thick rutile-like 
nanowires. 

In order to check the reliability of the DFT-DOS profiles and transport gaps, we have carried 
out more accurate electronic structure calculations for the rutile-based nanowires with 4, 12 and 
16 functional units. The choice of these testing structures is based on: i) their transport gaps are 
easy to identify (see [6]), while the thicker nanowires show states inside the gap region; ii) their 
small cell sizes lower the substantial computational effort of such calculations. For this purpose, 
we have calculated the many -body corrections using two different formalisms: a) quasi-particle 
corrections have been calculated within the GoWo formalism^ 41 * 42 * (quasi-particle Green function 
screened interaction approach) to properly account for the exchange-correlation self-energy; b) the 
GW + Bethe-Salpeter equation (BSE) approximation,^^ where the quasi-particle GW approach 
is further corrected to account for excitonic effects (electron-hole interation) in the unoccupied 
and excited states. The GW + BSE formalism includes, by construction, excitonic effects, which 
are important in highly correlated systems, and combines the excellent performance of the GW 
+ standard DFT approach for the ground state electronic structure with the high accuracy of the 
BSE methodology for the calculation of excitations. These calculations have been carried out 
with the Yambo package. ESI More than 200 unoccupied electronic bands were necessary to obtain 
converged self-consistent results. 
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[2] shows the transport gaps (in eV) predicted by the Go Wo and GW + BSE formalisms for the 
three rutile-like nanowires. The comparison between the DFT and Go Wo gaps indicates that the 
first-order quasi-particle corrections are small: the Go Wo gaps are larger than the DFT gaps by 
0.27, 0.14 and 0.1 1 eV for the nanowires with 4, 12 and 16 Si02 units, respectively. The excitonic 
corrections are larger. These lead to an additional increase of the gap by 1.18, 0.60 and 0.28 eV 
for the same three nanowires, respectively. The percentage change of the gap due to self-energy 
corrections, 

gG W _gDFT 

AEf FT -* G ° w °(%) = 100 x 8 8 , (6) 

is nearly independent of wire diameter (about 6 %). On the other hand, the additional percentage 
change induced by excitonic effects, 

F GW+BSE _ F G W 

AE G W^GW + BSE {%) = 1(X) *g _Eg (7) 

8 v I pG W > v ' 

decreases as the nanowire diameter increases: 25 %, 20 % and 16 %, for the three nanowires, 
respectively. Thus, the largest gap corrections due to many-body effects occur for the thinnest 
nanowires. It is worth noticing, however, that the trend of a decreasing gap with increasing 
nanowire diameter predicted by the standard DFT calculations is maintained after including many- 
body corrections. 

Electron distribution and charge transfer 



In order to understand the bonding mechanism we have carried out an analysis of the electronic 
charges of the atoms. Theoretical calculations giving the electronic charge distribution in the 
system can, in principle, provide this type of information, but the way to extract it is not unique. 
The usual output of quantum mechanical calculations is a continuous electronic charge density 
and it is not evident how to partition electrons amongst fragments of the system such as atoms or 
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Table 3: Average s and p orbital charges (in e~ units) in the Si and atoms, and net electronic 
charge gained / lost (in e~ units), for the /3-cristobalite, ortho-tridymite and the rutile-like nanowire 
with 24 functional units. Notice that the values shown correspond to one Si atom and two O atoms, 
that is, to a SiC>2 unit. 





/3-cristobalite 


Ortho-tridymite 


Rutile 


Si atom 


s orbital charge 


1.1 


1.0 


1.0 


p orbital charge 


2.2 


2.2 


2.2 


Net charge gain / lost 


-0.7 


-0.8 


-0.8 


O atoms 


s orbital charge 


3.4 


3.4 


3.3 


p orbital charge 


9.3 


9.4 


9.5 


Net charge gain / lost 


0.7 


0.8 


0.8 



molecules. Different schemes have been proposed, some of them based on localized electronic 
orbitalsP^EEl and others based on the electron density.^ 

Taking advantage of the localized basis set used in the Fireball code, we can evaluate the 
electronic charge in each localized orbital at the end of the self-consistent calculation.! 19 ! 22 ! 23 ! The 
Lowdin method^^has been employed to obtain the orbital charges. With those values, the sign 
and the magnitude of the net charge transfer between atoms can be obtained. In order to check 
the accuracy of the results we have compared the atomic charges with the Bader atomic charges,^ 
obtained using electron densities from the plane- wave DACAPO code. There is excellent agreement 
between atomic charges calculated in the two ways. 

Electronic charge is transferred from Si to O atoms due to the higher electronegativity of oxy- 
gen. [9] shows the average net charge transfer, in e~ units, from Si to O atoms, per SiC>2 unit. In 
other words, this is the average charge transfer per Si atom. Charge transfers in the rutile-like 
nanowires range between 0.73 and 0.81 e~ per Si atom. These are slightly larger than the charge 
transfers in the cristobalite- and tridymite-like nanowires, 0.68-0.73 e~. The larger charge transfer 
in the rutile-like nanowires correlates with their higher stability. For this purpose, one can notice 
the qualitative similarity of the data in [4] and [9} 

[3] shows average s and p orbital charges, in e~ units, in the Si and O atoms for three nanowires: 
/3-cristobalite and ortho-tridymite wires, and the rutile-like nanowire with 24 Si02 units per cell. 
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Figure 10: (Color online) Bars diagram showing the s and p orbital electronic rearrangements in 
the Si and O atoms of the /3-cristobalite-like nanowire (in e units), per Si02 functional unit. The 
net charge transfer is also shown. Blue and red filling indicates electronic charge lost and gained, 
respectively. 

The net electronic charges lost (by each Si atom) or gained (by two O atoms) are also given. The 
results for the two tridymite-like nanowires are very similar. Also the results for the six rutile-like 
nanowires are very similar. This, plus the information in|3j indicates that the average atomic charge 
populations, total and orbital-resolved, are similar in all the nanowires. We then conclude that the 
chemical bonding in these nanowires exhibits common features, quite independent of the detailed 
structure of the nanowires. 

The chemical bonding can be described as a mixture of covalent and ionic contributions. In all 
cases, the p-like electron population of the O atoms gains 0.70-0.75 e per O atom (1.4-1.5 e~ per 
Si02 unit). The Si atoms also increase their p-\ike charge by 0.1-0.2 e~ per atom. This general 
increase of p-like population in the nanowires arises from a loss of charge in the s-like orbitals 
(0.9 e lost by the Si atoms and 0.30-0.35 e lost by each O atom). This is summarized in the bar 



diagram of 10 for the /3-cristobalite nanowire. 



In order to visualize the spatial distribution of the electronic charge, 11 shows surfaces of 



constant electron density in 3 dimensional perspective for the rutile-like nanowires with 4 and 12 

_ ° -3 

Si02 units per cell. The densities selected are 0.1, and 0.03 e A . The results are similar for 
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4 functional units 12 functional units 

TOTAL DENSITY 

0.1 e/A 3 0.1 e/A 3 




0.03 e /A 3 0.03 e /A 3 




Figure 1 1 : (Color online) Surfaces of equal electron density for two rutile-like nanowires with four 
and twelve Si02 units per unit cell. The densities selected have the values 0.1 and 0.03 e A for 
the total densities, and 5- 10 4 e~A for the HOMO and LUMO densities. Oxygen and silicon 
atoms are represented by red and yellow spheres, respectively, and the iso-density surfaces by grey, 
blue and green shiny tones for the total, HOMO and LUMO densities, respectively. 
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the two nanowires (and for others not shown here). The surfaces corresponding to a density of 
0. 1 e A are formed by a set of disconnected pieces showing a concentration of charge around 
the oxygen atoms, aided by the charge transfer from Si to O atoms. The surface corresponding to 

_ 9 -3 

the lower density, 0.03 e A , is a connected surface displaying the mixture of ionic and covalent 



bonding in the nanowires. 1 1 also shows the charge densities of the HOMO and LUMO orbitals, 
again for the rutile-like nanowires with 4 and 12 Si02 units per cell (all with values of 5-10 -4 

_ ° -3 

e A ). In the thin nanowire, the HOMO orbital is localized near the planes formed by the most 
external oxygen atoms (each of those planes contains four oxygen atoms). In contrast, the LUMO 
orbital is localized near the planes containing the distorted squares formed by four Si and four O 
atoms. For the thicker nanowire, the HOMO orbital is placed in between the planes containing the 
Si atoms, and the LUMO orbital mostly surrounds internal Si atoms. 



SUMMARY 

This work presents ab-initio simulations for Si02 nanowires. The structures of the nanowires 
have been determined: i) by relaxing nanowires cut from Si02 bulk crystals (/3-cristobalite, ortho- 
tridymite and /3-tridymite), and ii) by relaxing nanowires with a rutile-like structure. The analysis 
of the electronic states indicates that many of these nanowires have a semiconducting character, 
different from the insulator character of the macroscopic bulk-crystals. However, for the thick 
rutile-like nanowires we found the gap congested by the presence of surface states. The bonding in 
all the nanowires can be described as partly ionic and partly covalent. Novel uses for nanoscopic 
Si02 nanowires may arise from insights gained about the electronic density of states, chemical 
bonding, charge density distribution and cohesive energies per Si02 unit as a function of nanowire 
thickness and structure. 
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